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The field-theory for multifractals in percolation is reformulated in such a way that multifractal 



exponents clearly appear as eigenvalues of a second renormalization group. The first renormalization 
group describes geometrical properties of percolation clusters, while the second-one describes elec- 
trical properties, including noise cumulants. In this context, multifractal exponents are associated 
Oh. with symmetry-breaking fields in replica space. This provides an explanation for their observability. 

It is suggested that multifractal exponents are "dominant" instead of "relevant" since there exists 
an arbitrary scale factor which can change their sign from positive to negative without changing the 
£SJ ■ Physics of the problem. 
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ON \ I. INTRODUCTION 

O 
ON 

The renormalization group and critical phenomena have provided, over the years, the key concepts which allow us 
' to understand problems ranging from phase transitions to percolation. Even though infinite sets of exponents, such as 
ON ■ crossover exponents, were calculated early after the introduction of the renormalization group in critical phenomena, 
jjj attention is usually focused on a few relevant exponents. This focus on a few exponents is justified since observable 
quantities in general couple to many renormalization group eigen-operators, including the most relevant ones, which 
eventually dominate their behavior. That lore was challenged relatively recently by the appearance of infinite sets of 
measurable exponents in various fields. Problems where such infinite sets occur are collectively known as multifractal 
problems, § || even though some of these problems have a quite different physical nature. We discuss here the 
q \ problem of electrical properties of percolating networks. || (7| 

The difficulties associated with the formulation of a Lagrangian field theory for multifractals have been first discussed 
by Ludwig and Duplantier. j|] 0. In particular, the exponents x n describing multifractal behavior are a convex 
function of n, while the analogous exponents of field powers in a field theory are in general concave, as a consequence 
of stability and correlation inequalities. ||] But, as stated in Ref. ||, powers of composite operators, like derivatives, 
can exhibit multifractal behavior. It is therefore possible to formulate for multifractals a field theory which is amenable 
to a renormalization group analysis. This has been done for percolation by Parks, Harris and Lubensky JTc| (PHL). 
This field theory however does have some peculiarities which make it different from usual field theories. In the present 
paper, we reinterpret the field theory of PHL so as to make special features clearer. In fact there are also some 
differences in the way we set up and interpret the field theory. A detailed discussion of the differences and of the 
reasons which motivate our approach have appeared. fll|| The structure which will emerge here is that geometrical 
properties of percolation clusters are described by a standard field theory on which one can perform a renormalization 
group analysis. We call this the "first renormalization group". By contrast, multifractal properties, originating from 



the electrical properties of the cluster, are described by a "second renormalization group". [[L2| |13| The structure of 
this second renormalization group depends on the first renormalization group: Once the usual recursion relations for 
percolation are derived, the recursion relation for each multifractal moment must be found by a further projection onto 
an appropriate eigenbasis. In the second renormalization group, the role of the fields is played by the microscopic noise 
cumulants v s . These fields v s are conjugate to replica-space gradients of the operators, for example k 2s <I>k(q) < E > -k(— q) 
where q is defined in the usual Brillouin zone while k is the replica-space extension of the Fourier variable for the 
voltage. Without changing the Physics, the scale of k can be changed by an arbitrary factor. The existence of this 
arbitrary scale factor is a direct consequence of the linearity of Kirchhoff 's laws. This has no equivalent in usual 
critical phenomena. Due to this arbitrary scale, the multifractal exponents can change sign, rendering inappropriate 
the usual classification as relevant or irrelevant. We suggest, therefore, to call the multifractal exponents "dominant" 
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exponents since they determine the leading scaling behavior of observable quantities while corrections would come 
from "sub-dominant" exponents. The fields for multifractal moments are associated with operators which break 
rotational (permutation q|) symmetry in replica space. This seems to be the reason why the multifractal moments 
are each associated with a different dominant exponent, and not just to a few ones. 

In the following section, we make a short review of the phenomenology of multifractals in percolation. Then we 
derive the field theory and proceed with the renormalization groups. The discussion section clarifies some points of 
interpretation. 



II. PHENOMENOLOGY OF MULTIFRACTALS IN PERCOLATION 



For completeness, let us recall how the infinite set of exponents appears in percolation through the problem of 
noise. Suppose that the conducting resistors of a percolating network are fluctuating independently in time. The 
total resistance of a given network is then a random variable in time whose cumulants depend on those of each 
component resistor. The cumulants of a given order are assumed to be the same for all component resistors. The 
cumulants of the total resistance R are, in principle, accessible experimentally, and measurements of the second 
cumulant, corresponding to 1// noise, have actually been performed. [Q Theoretically, for systems of finite size L at 
bulk criticality, one finds, after averaging over the microscopic noise, that the cumulants, C|j (L), of order n scale as 

< C%> (L) > c = v n < J2 >c~ L~*« (1) 

a 

where C represents average over percolating lattice configurations, i a is the current that flows in branch a of the time 
averaged network and where v n is the amplitude of the n'th cumulant of the elementary resistance fluctuations. For 
example, the usual electrical noise amplitude for one microscopic resistance r is obtained from 

v 2 = {SrSr} f (2) 
where {}^ refers to time average over the noise. The first equality in Eq. ([!]) follows from Cohn's theorem 

SR = J2 Sr <* l l ( 3 ) 

a 

where the total input current, Ii n j., is unity. Each exponent x n is different and is not a simple linear function of n, as 
commonly occurs in critical phenomena under the name "gap scaling" . JjJ Such an infinite set of exponents also arises 
in the other analogous problems mentioned above. In all these problems, one is assigning to parts of a fractal network 
a weight, or a measure, which is obtained from the solution of Laplace's equation. Here that weight corresponds to 
the power dissipated in a bond. 

Finally, let us recall that the positive integer moments, < ^ Q i 2 ™ >c, suffice to characterize completely the 
distribution of the currents flowing through the network. |p.5| But, as stated previously, the exact value of the integer 
moments is necessary to reconstruct all the information so that the leading scaling behavior of the positive moments 
does not suffice to find, for example, negative moments. In the following, we concentrate only on the scaling behavior 
of the positive integer moments. Note that in order to keep the same notation as PHL JlCfl , we work with exponents 
defined by 

tp s /u=-x n (s = n). (4) 



III. FIELD THEORY FOR THE GENERATING FUNCTION 

A. Generating function for the cumulants 

Since Kirchhoff 's law can be obtained by minimizing the entropy production, or the dissipated power in the electrical 
network, it is natural to start from a generating function for the resistance between two points y and y' of the network 

W (k; y, 2/0 = / d [V(x)] exp (-H + ik [V(y) - V(y')}) (5) 
-Mv min 
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where AVmin and AV max are, respectively, the minimum and maximum voltage drops for a finite size system, and 
where H is given by, 

H=\ <r b (x,x';CJ)[V(x)-V(x')) 2 (6) 

<x,x'> 

with <7f,(x, x'; C, /) the conductivity of the bonds a linking each pair of points x and x' for a given configuration C of 
the random resistor network. For each given configuration of the random network, the conductivities fluctuate in 
time. In other words, 

<ib(x, x'\ C, /) = <jq{x, x'; C)(l + e), (7) 

where e < 1 is a random variable whose probability distribution is given by / (e). Hence, there are two types 
of averages to perform: The usual bond-disorder average and, for each lattice configuration, an average over the 

microscopic noise. To obtain a Hamiltonian with the same structure as for a spin system, (e.g. J (x — x S(x)S(x )) 

one uses Fourier transforms. Since the potential differences are on a bounded interval, determined by the boundary 
conditions, it is possible to use discrete Fourier series. Formally, we may write 

H= Vb(x,x';CJ)J2Ake lk(v ^- vix '» (8) 

<x,x'> k 

where the Fourier coefficients are given by 

1 r AV max 

^ k = "atT / d ( AV ) [AVfe lfcAV \ (9) 



AV m 



Because the k = terms correspond to a uniform distribution of voltage, they will be discarded. In the theory of 
Stephen, [fl6[ the order parameter tp is the defined in such a way that its autocorrelation function vanishes in the 
non-percolating phase and decays exponentially in the percolating phase. This order parameter is 

My) = e* kV{y) (io) 

where the wave vector k takes a discrete number of values 27rn/(AV max ). The presence of a lattice, instead of 
a continuum, leads to the existence of an ultraviolet cutoff corresponding to the minimum voltage [[lTj = 
2ir / (AV m - m ) . That cutoff becomes infinite in the percolation limit er^ 1 — > 0. For a given realization, we may then 
write H in the form advertized, namely 

H =E A fe E <Tb{x,x?;CJ)il> k (x)MJ) (11) 

k <X,x'> 

where, in the second sum, we consider only the connected links of the network. To understand what follows, it is 
important to realize that the fc's scale as AV max _1 and that we must use a discrete Fourier expansion, since the 
replica method is only valid for finite system. The infinite-system limit is taken at the very end. 

In the limit where u^ 1 — > 0, the saddle point configuration for the voltage drops obeys Kirchhoff's laws and, in this 
limit, we obtain for the generating function 

W(k;y,y')=<My)^-k(y')>= Z e -i^»'^ (12) 

where 

r -k J2 °b(x,x';Cj) [v(x)-v(x')} 2 

Z = d[V(x)} e <*■*'> = / d[V(x)]exp(-H) (13) 



and R(y,y';C,f) is the resistance between nodes y and y' for a given configuration C of the random network, and 
realization / of the noise. As in PHL 0] the cumulant averages for the resistance noise may be obtained from the 
generating function. We first average Eq. (^|) over the probability distribution for the noise / and then we average 
over the disorder as follows, 

/jexp(-ifc 2 i?(y,y';C)))j^ v{y,y';C)\ (14) 
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where the function v plays the role of a conditional probability that the two points y and y' are connected. Expanding 
the left-hand side in cumulants of the resistance, Eq. (O) becomes, 



exp 



E 

S>1 



2 s s\ 



(15) 



c 



The cumulants of the resistance fluctuations C^\y,y'; C) may thus be obtained from derivatives of the generating 



function since, by Eq.(Q) and in the a Q 1 — > limit, they are proportional to the microscopic cumulants, fa] 

(y,y';C)) ! ,(y, y ';C')) c 
»[({<*k(»)*-*(» , ))} / > 



(-l) s 2 s s!fc- 2s -. 



5 V s 



(16) 



l«i=o , v; 



Thus, the macroscopic cumulants of the noisy resistor network can be obtained from the autocorrelation of the order 
parameter, after averaging over bond disorder. It is important to note that the quantities v s play the role of fields 
whose conjugate operators will contain polynomials in fc, as explained later. 



B. Replica method and effective action 

Since there are two types of averages over random variables, we introduce two types of replica, as suggested by 
PHL |[o]]: N replica for the average over noise (/), and M replicas for the average over bond disorder (C): 

n M 

{z N } f } 

^ Z h) c = V " Nm' C 7 (17) 

As usual, the limiting process is justified only for finite systems. In other words, the limit L — ► oo is taken after the 
limits N, M -> 0. 

We do not repeat the details of the derivation of the field theory. [|l0| As usual, it proceeds by introducing Hubbard- 
Stratonovich variables <& conjugate to each i/j appearing in Z NM . Expanding in powers of tp and performing the 
integrals over the voltages appearing in the original generating function, one generates a power series in $ which 
can be re-exponentiated to yield an effective action for the <&. These variables are now the operators of the field 
theory. For Hubbard-Stratonovich transformations, the generating function < ip-\ SL (x) , ip--k(x) > is simply proportional 
to < $k (x) ^-k (x) >, hence all we need is the effective action for the operators 3> k 



£($) = !/ d d x Afc £ k r k [$ k (x) $_ k (x) + V$ k (x) • V$_ k (x) 
+ f / d d x (Afc) 2 £ kl+k2 ^o * kl (x) $ k2 (x) $_ kl _ k2 (x) . 



(18) 



The k in Eq. @ is any of the k Fourier variables whose components are labeled k a p in the NM— dimensional 
replica space and!? E^ o(? fc^ is the square-modulus of k. There are (L MN — 1) operators < I > k (a;) at spatial point 
x (the k = case is omitted). The k a p are conjugate to the electrical potentials of the replicated systems. They 
contain a scale factor correponding to the scale of the electrical potential. In the limit of geometrical percolation 
(i.e. no transport property), the k a /3 are all zero. Note that despite the notation, the k are tensors of rank one and 
not two, as far as "rotations" in replica space are concerned. As usual, to find critical exponents the k dependence 
of the coupling constant M3 can be neglected, but that of r k is crucial. The scaling of the terms entering r k can be 
infered from the scaling properties of the generating function. Writing explicitly the dependence on the microscopic 
cumulants v s , one obtains for the scaling properties of the generating function, 

G k (y - y',p-Pc,{ v s} s >i = G Ak (y-y',p-Pe, { A ~ 2s M s >i ( 19 ) 

since k scales as -yrr — ■ Indeed, the cumulants C^'(y,y'; C) are linearly proportional to the microscopic cumulants 
v s so that, as can be seen from Eq. (|l5|), the autocorrelation function for the order parameter does have the scaling 
property (|19|). Expanding for small v s , the most general form for the r k is then an homogeneous polynomial in k and 

v.. 



4 



s>1 



7-k =P-Pc + ^k^ s ( ^ Uai«s 3 ...U Sj P( Sl .s 2 ,...,s z ) (-Ac, >(3 , ...) J . (20) 

S2 + --- + Si=S / 



In Eq.(|20|) the functions P s (..., 6*k, ■■■) depend on the angular variables Ok in the n = MN dimensional replica space. 
For each s =(s\, Sa, •••) Si), these polynomials can be expressed as a linear combination of the spherical harmonics for 
the MN angular variables 0^ . Their explicit expressions depends on the particular distribution / for the noise of 
the elementary bond resistor. The expansion in Eq. (^0|) is justified by the fact that for a finite system, the effective 
action is an analytic function of k and hence, by scaling arguments, of v s . For s = 1, it will be seen in the Appendix 
that P\ (#k) = 1, i.e. k 2 is the only polynomial of degree 2 which contributes and it is an eigenpolynomial of the 
second renormalization group. 

Finally, observable quantities arc obtained using the standard replica-method identity, 

G k [y- y',P -Pc, {vs} s >i > J mi-) = ({(Tpk(y)ip-k(y')}} f ) c , 21 n 

- limjv^oA/^o ($ k (y)*-k(y')> 

In the absence of symmetry breaking in replica space, the modulus of the replicated k in Eq. ( pl| ) is equal to fc in 
expression (|l5|). 



C. Scaling properties 

Near the percolation critical point, the generating function scales as, 



lim r -i „ 

L — >oo ; cr n — >u ; n- 



+0 G k (y ~ y' ',P - Pc,{Vs} S >1 .Anj.) 



lim L^oo ; .-^0 ; n^O A 2 ~' 7 »' _<i Gfc ((j, - „') A -1 , (p - Pc ) A 1 /^ , { Vs \i>./*, } ^ , 



(22) 



which also implies 



lim r „ 



+ ; rwO °k (V -V',P- Pc RL>1 i knj) = 

(p - p c f S ((y - t/') (p - Pc r , {«, (P - Pc)-^} s > x ) 



(23) 



where S 1 is a scaling function. As usual, the correlation length £ behaves as (p — p c ) . The L — > oo must be taken at 
the end for the replica method to be valid. From the way the problem is set up, the n — > limit must be taken before 
the Uq 1 — > limit. However, as we shall see in the following sections, we need to add another condition to linearize 
the second renormalization group equations, namely that 

v s « v\ ; Vs > 1 

This inequality is consistent with the phenomenology of 1/ f noise. 



IV. RENORMALIZATION GROUP APPROACH: 



A. First, and second renormalization group. 

The renormalization group equations are obtained JL0[ | from the usual procedure for a cubic Landau-Ginzburg 
functional. Using the Wilson approach, operators ^(q) whose wavevector q is in a shell A/6 < q < A, near the 
cut-off A coming from the physical lattice, are traced over. It is always possible to choose the lattice spacing units 
such that A = 1. Lengths are then rescaled by b, while the operators are rescaled as follows: 

*k(q' = M = & (_<i ~ 2+? ' p)/2 $k(q) (24) 

Let Kd be the surface of the d— dimensional sphere. To one-loop order, in dimension d = 6 — e near the upper-critical 
dimension d = 6, the differential recursion relations for rk and for the coupling constant g = Kdu\/2 take the form 
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^= (2-77 p )r k - ff E k (25) 



= (e - 3?7p) g 



where I is defined by 6 = e l . In the first equation, £ k is the self-energy-correction to one loop |lC| ] 

S k = lim n ^o(Afe) Ep, P +kjio 

G [p] G [p + k] . 

= -2G [k] G [0] + E k 1 j 

Since we are looking for a linear renormalization group in v s , it suffices to restrict ourselves to the Green's functions 
expanded to first order in v s . Furthermore, the momentum-shell integral is for q 2 = A 2 = 1 so that the Green's 
functions appearing in ( p6|) are of the form, 

Glp + kf 1 = l+p- Pc + Y / Vsk 2s P$(...,9 k ,...). (27) 

S>1 

The usual percolation fixed-point, v s = 0, and p = p c describes the geometrical properties of percolation clusters. 
These properties are completely independent of the electrical transport properties of the network. In other words the 
field tq is the coefficient of a polynomial of order in k so that the calculations of the fixed point and of the usual 
geometrical exponents r] p v p are independent from the values of the fields v s which are all associated with higher-order 



polynomials in k times the field operator $ k $_ k . Hence one obtains the same results as in ref. 10 namely g* = e/7 



for the fixed-point value of the coupling constant and rj p = — e/21, ^p — i + ff f° r the exponents. These exponents, 
and fixed-point values of g and r$ come from what we call the first renormalization group, whose predictions relate 
to the geometrical properties of the percolation network, and hence are independent of electrical properties. 

Going further to obtain the exponents ip s associated with v s , one has to linearize Eq. (p5| ) in v s . To the order 
we are working, the v s are all multiplied by the same power $ k , namely "J 2 ., but by different powers of k, namely a 
homogeneous polynomial of order 2s. Hence, we consider the renormalization group as a renormalization equation in 
the space of homogeneous polynomials in k. This is what we call the " second renormalization group" . 

To be more specific, let us start with the simple case of d larger than the upper critical dimension (d c = 6). The 
renormalization group equations for r k then read, 

(28) 

In other words, the recursion relation for the v s is obtained by linearizing around the fixed-point v s = 0, s = 0, 1... 
To first order in v s , one has 

r k =p-;p c + ^k 2 V s (...,0 kQ/J ,...) (29) 

S>1 

and thus the scaling in Eq.(p4|) leads to, 

dv„ 

-^ = 2,, (30) 
Above six dimensions then, any P(k) = k 2s P s (...,# k ...) is an eigenvector of the second renormalization group. 



The eigenvalues are all identical, as is already known. |5j Thus, even though the precise form of -P s (k) depends on the 
microscopic noise distribution, the eigenvalue for v s does not. Below six dimensions, things are less trivial since the 
self-energy E k depends functionaly on P s (k). Nevertheless, we will show that there exists an eigenbasis of polynomials 
in k whose eigenvalues are labeled only by the order of the polynomial. This means that whatever the starting P s (k), 
the eigenvalue for v s depends only on s. 

Below six dimensions, the RG equations must be obtained by linearizing the self-energy S k as a function of v s , s > 2. 
We show in the Appendix that the eigenpolynomial basis is obtained, in the limit v s << Vi, by solving for each integer 
s the following equation 



-2 r+°° du f +o ° dt exp {- (u + t)(l + r )} ' 



t+u 



{exp(^- 2 A)p(k)}=4P(k) 



(31) 



which correponds to the limits a^ 1 — > and n — > for the self-energy. In Eq. ( |3l|) the Laplacian operator A acts 
on k and its limit, as n — > 0, is discussed in the next section. Because r k in Eq. (Eq) is linear in v s , the eigenvalues 
for the "second renormalization group" are obtained as 2 — r] p — g* (2 + d s ). Everywhere below, g will be taken at its 
fixed-point value, g* . 



G 



B. Eigenvalues and eigenvectors of the second renormalization group. 



Eq. (|25|) and (|31j) define the renormalization group equation for homogeneous polynomials of order 2s. To diago- 
nalizc this equation we need cigcnfunctions for the Lapalacian operator. In n-dimensions, the most general form for 
these eigenfunctions is given by pljfl 



k~ — J, 



2r+- 



ilk) Y 2r (...,6 k) ...) 



(32) 



where Y 2r (...,6i, ...) is a spherical harmonic of order 2r in n = NM dimensions, J r is the Bessel function of order r 
with eigenvalue — 7 2 . To expand polynomials of order 2s on this basis, we first recall that the maximum value of the 
replica-space k = VI? is given by the ultraviolet cutoff Ak = t max = 27r/AV m ; n = ao2ir/i m i n . We can then proceed 
as follows. First note that an arbitrary polynomial of order 2s can be written as, 



P s {k a0 ) 



0<r<s 



Q>s,r Y 2r (• 



•) 



(33) 



where, to simplify the notation, we do not write internal indices related to the degeneracy of the spherical harmonics. 
Then, we write 



k 2s Y 2r (..., k , ...) = (k 2s - k 2r A 2(s - r) ) Y 2r (..., k , ...) 



\t 2r Af s - r) Y 2r i 



(34) 



The last term is an eigenfunction of the Laplacian with eigenvalue zero, while the first term can be expanded using 
the Fourier-Bessel expansion, |2C|] which is uniformly valid for functions with zero value at the end of the interval. 



(' 



-2.s 



2r A 2(s-t 



k ar A 



Y 2r (. 



Y,i>ihk~^J 2r+ ^ (^0 )Y 2r { 



■) 



(35) 



where {Ct} l>1 are the zeros of the Bessel function of order 2r + (n — 2)/2. 

We can now substitute this expansion in the kernel of the integral appearing in the recursion relation Eq. (^ 
The contribution of a general term of the series will have the form, 



du 



dt exp 



2v\ 



0_ 
A fe 



J. 



2r+'- 



A fc 



Y 2r (. 



•) 



(36) 



The n = NM — > limit is obvious. Taking the a 1 — > limit simplifies further the equation considerably since, in 

the argument of the exponential, [«iA|] 1 is proportional to a^ 1 . This means that effectively every term of the series 
behaves as if it had the same eigenvalue for the Laplacian, namely 0. Since k/A^ is bounded between and 1, the 
series can be resummed, and any polynomial of order 2s is an eigenvector with eigenvalue [fiof 



d s = -2 



+ OG 



du 



+ OC 



dt exp (— (u + t)) 



t 



-i 2s 



t + U 



(2a+l)(a+l) 



(37) 



The main conclusion of the section is that the space of homogeneous functions of order 2s can be characterized by a 
single eigenvalue which depends only on the order of the polynomial. The form of the microscopic distribution of the 
noise is, therefore, not relevant for scaling properties of the macroscopic resistance fluctuations. Only the condition 
v s << vx , Vs > 1 which was used to derive Eq.(31) needs to be satisfied. The critical exponents associated with the 
fields v s are thus 



tp s = u p (2 - r] p - g* (2 + d s )) = 1 

I Lubensky [y v = \ 
applies also for the case s = 1, as written above. 



1,2.. 



(38) 



as found by Parks, Harris and Lubensky (y v = ^ + ; r/ p 



7(s + l)(2s + l)' 
= — ^j-). Jl0[ It is shown in the Appendix that this formula 
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V. DISCUSSION 



A. Gap scaling 

The scaling of the usual thermodynamic observables is normally trivially obtained from a few exponents only. This 
is usually refered to as "gap scaling". Gap scaling also occurs in the present case. For example, 

^l^j^o-k-.'r/- (39) 

That this applies to multifractals in percolation was verified by numerical simulations in EL Clearly, one can also 
define universal amplitude ratios, jla] pl| 



B. Symmetry breaking 

With v s — for all s, the action is invariant under the global transformation 

$k = (40) 

where 1Z is a rotation (permutation) of the vector k in the replica space of dimension MN. In other words, the action 
transforms according to the unit representation of the group 0(MN). When v s ^ 0, that symmetry is broken, since 
polynomials of higher degree transform like higher-dimensional representations of the group O(MN). |2^] To have 
every v s associated with a different representation of the symmetry group is a necessary but not sufficient condition 
to have an infinity of observable exponents. Indeed, operators of different symmetry could couple when higher-order 
corrections to the e expansion are evaluated. All this is analogous to what happens with symmetry-breaking fields in 
other critical phenomena models, such as the XY model for example. ^3|. The exponents tps here are analogous to 
the crossover exponents ip„ of the XY model. p4j. 



C. Dominant exponents 

At first sight, the perturbations associated with the v s are all relevant since all the exponents tps are found to 
be larger than zero. There are two important differences however with critical phenomena (say the case of the XY 
model) 

a) There is no physical realization that we know of for the lower-symmetry fixed-point towards which the system 
rescales when one of the symmetry-breaking perturbations is different from zero. All physical observables are deriva- 
tives evaluated at a zero value of the symmetry-breaking fields v s : In other words, the exponents ?p s are crossover 
exponents associated with the symmetric fixed-point. 

b) There is an additional freedom to rescale k at each iteration as seen in Eq. (|l9[) . This allows one to formulate 
the renormalization group in such a way that only a finite number of operators are relevant! Indeed, for the usual 
percolation fixed-point, the rescaling of the operators is found by choosing that the coefficient of the spatial gradient 
term in Eq. ( |l8| ) to be a constant. Since the recursion relations for U3 and tq are completely independent of k, the 
geometric percolation fixed-point is the usual one. The scale factor for k, by contrast, may be chosen at will. This 
influences the recursion relations for the v s and hence the corresponding %j) s exponents. More specifically, the rescaling 
part of the renormalization group transformation may be written as follows 

*k (q) — $U«k W = qb) = & (-<i-2+^)/2$ k {q) . (41) 

As an example, we choose the scale factor a such that the total resistance is kept constant under rescaling of all 
lengths by a factor b. This is done by first noting that after eliminating the degrees of freedom, the scaling of v s is 
obtained by keeping the corresponding terms of the action invariant 

v y'(2+d a ) k ,2s {ql)] 2 (Afc ^MiV ddql = Vsk 2s [$k {q)] 2 {Ak) MN ^ (42) 

In the limit NM — > 0, (Afc)^ 7 ^ does not come in since the infinite-system limit (Ak = 0) is taken last. Setting a = b v , 
and b — e l , we obtain, 
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-f = [2- Vp -g* (2 + d s ) + 2sy)v s . (43) 



Choosing 

y = -\[2-Vp-9*{2 + d x )] (44) 

to keep the voltage across the network constant (i.e. v\ does not scale), the exponents associated with v s become 

ip a -> ip s - 2s -01 (45) 

so that the fields v s now appear irrelevant for s > 2. In other words, since the ip s are a decreasing function of s, we may 
always choose the scaling dimension of k such that only a few of the ip s exponents are positive, without influencing 
the Physics. These statements can be rephrased in a more physical way by recognizing that the field theory of PHL 
corresponds to computing the power dissipated between points x and x' when a unit current is injected between these 
points, whatever the distance between x and x': One could just as well decide to rescale at unit voltage instead of 
unit current, and this would correspond to multiplying k by a scale factor at each iteration. The size dependence of 
Afe does not matter, since, in the limit o~ x — > 0, A& goes to infinity independently of the system size. The rescaling in 
k is associated with the scaling of applied voltage so that the scaling in k and q space are independent of each other. 

The remark of the last paragraph may also be formulated as follows: The analysis that we have done to find the 
eigenpolynomials for Eq.(^5|) shows that the latter equation is like a "second renormalization group" with fields v s 
which describe the electrical properties of an object whose (critical) geometrical properties are given by the "first 
renormalization group" , with fields rg and U3 . The first renormalization group has properties totally independent 
from those of the second while the second is slaved to the first. The v s are conjugate to replica-space gradients of the 
field operators $. In other words, they are conjugate to polynomials of degree 2s in k times The rescaling of k 
in the second renormalization group is arbitrary, and this is fundamentally due to the linearity of Kirchhoff 's laws. 
Instead of referring to relevant or irrelevant exponents, for that second renormalization group it makes more sense to 
call them "dominant" exponents since one expects that the observables which are connected to ip s are also coupled 
to other operators giving corrections to scaling (sub-dominant exponents) . 

D. Observability 

The ip n of the XY model are not all relevant exponents. In fact, for n > 4, they correspond to irrelevant operators. 
|jj [^3| [^4) (2^] While they are only a subset of all possible irrelevant exponents, they are, however, special because, 
for increasing values of n, they represent the leading scaling behavior of operators with lower and lower symmetry. 
It is their symmetry instead of their relevance which seems fundamental for their observability. The same remark 
applies for the "dominant" exponents tps/v = —x n discussed above. 

VI. CONCLUSION 

We have shown that the field theory for multifractals in percolation has a special structure which allows multifractal 
exponents to have properties which do not usually appear in standard critical phenomena. They follow from symmetry- 
breaking operators in a "second renormalization group" to which an additional normalization freedom (e.g. scaling at 
constant voltage or constant current) is associated. This freedom allows one to arbitrarily shift the crossover exponents 
(while maintaining the observable quantities unchanged). We propose to call these exponents dominant since, even 
though their value can be shifted, they are trivially related to the leading scaling behavior of operators characterized 
by a given symmetry. A similar structure with a "second renormalization group" also occurs in dynamical systems 
p2| , and probably also in the field of localization. 
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APPENDIX A: RECURSION RELATION TO ONE-LOOP 



In this Appendix we derive the R.G. Eq. (|3l|). Let us recall that, to one-loop order, the recursion trelation for 
takes the form 



^= (2- Vp )r k -<?(2 + £ k ) 



where the self-energy is defined by 



Ek = lim 



+00 



— 00 j —00 



+00 



d Paf3 [G(p + k)G(p)] 



with 



[G (k)f 1 = 1 + p - Pc + £ Vsk 2s P§ (• • ., e k , ■ ■ •) 



(Al) 



(A2) 



(A3) 



Eq. (Al) is a non- linear recursion relation for the v s . We now show how one can linearize this recursion relation in 
the limit where n — > and — > 0. Using the Schwinger representation for the propagator 



+ 00 



G(k) = / du exp -u(G(k)) _1 



we get for Eq. (A2) 



lim 

n— >0 



p+oo p+oo p+oo p+oo 

/ du dt ••■/ da; Q/ 3exp{-u[G(x)] _1 -i[G(x + k)] _1 [ 

J JO J-oo J-oc ' 



By making the change of variables 



t 



t + u 



(A4) 



(A5) 



(A6) 



Eq. (A5) can be expressed as 



,. f + OO , f + OO , p + co p+oo , 

hm„^ Jo du J dt J_ oc ■■■ J_ ao dx 



a/3 



exp 



-uG 



tG 



(A7) 



For v s << vi, we can expand to first order in v s all terms in G which are not invariant under rotation (the only 
term which is invariant is k 2 and it is associated with the resistance). It is useful to define 



P^J>k 2s P ; 



(A8) 



s>2 



Therefore, Eq. (A7) reads as 



exp 



t r + oo , p+00 , p+00 p+00 , 

- lim^o Jo du J dt J_ oo • • • J^ dx af3 

- (u + 1) (1 +p - p c ) - vix 2 (u + 1) - fi^k : 



(A9) 



P 



t 

t+u 



p 



t 

t+u 



To interpret Eq. (A9) observe that in dimension 



dx i~J-™ dx n exp[ -vi(u + t)(xf + ... + 3%)]x£\..x£ i 

i vi(u+t) ) 



3 




1 


Pi- 


1 


2 p ipi< 


2i>i(«+i) 





(A10) 



10 



which gives asn->0 



n 



i 



l2P*pi\ I2vi(u + t) 



(All) 



pj P 



For analytic functions P, wc can write with k' = 

P(k' - x) + P(x + k') 

2 V V - r 2pi 

A l^i>o l^ Pl +...+ Pj =i {2p 1 )\...{2 Pi y.- h i ■■ 

The integrals can then be evaluated as 

lim^o d"xexp[-ui(w + i)x 2 ] [P(x - k') + P(x + k')] 



c 2pj JL 
J dk 



2P1 



(A12) 



^ J2l>0 {Avt(u+t)y Spj 
A 



r^2p 1 + ...+2p j - p 



2 [j^fej] ^ k ') = 2 exp [^t~ 2 A]P(k) 



(A13) 



where A is the Laplacian operator for the variable k. In this equation, we have used the fact that for s> 2 the term 
v\j^k 2 in Eq.(A9) cannot play any role to linear order in the v s . Now that we have interpreted Eq. (|A9|) , we can 
substitute it in Eqs.flA7|), flA5|) and flA^) to obtain, 



(A14) 



Iim o .-i_ >0 linin^o S k = 
lim^-i^o J2 s >2 2l >s Jo°° du Io° dtexp[- (u + t)(l + r )] 



u+t 



exp 



4di 



^(k) 



Because the thermodynamic limit is taken last in replica approaches, we can show, as is done in the text, that 
exp f^£~ 2 A has a well-defined limit as er^ 1 — >• 0. 

The case s=l can now be treated simply. Indeed, in this case, it is only the term vijj^k 2 which plays a role. 
Expanding it, the term linear in v\ in Ek is given by 



lim v\ I du dt exp f— (u + t)] k 2 

Jo Jo "t + u 



This shows that k 2 is an eigenpolynomial, as quoted in the text. Furthermore, the corresponding eigenvalue d\ does 
correspond to the s=l limit of Eq.(|38|). 



[1] M. Fisher, Universalities and Renormalisation Group theory, Lectures Notes in Physics , Spinger Verlag, Heidelberg, 186 
(1983). 

[2] U. Frish et G. Parisi, Turbulence and Predictability of Geophysical Flows and Climatics Dynamics. Proceeding of the 
International School of Physics Enrico Fermi, cours LXXXVIII, Varenna; T.C. Halsey, M.H. Jensen, L.P. Kadanoff, I. 
Procaccia and B.I. Shraiman, Phys. Rev. A 33, 1141 (1986). 

[3] C. Roland and M. Grant, Phys. Rev. Lett., 63, 551 (1989); P. Meakin, A. Conoglio, H.E. Stanley, T. Witten, Phys. Rev. 
A 34, 3325, (1986). 

[4] For a general review on Fractals and Multifractals see Fractals and Disordered Systems, Eds A. Bunde and S. Havlin 
(Springer- Verlag, New York, 1991), see also G. Paladin and A. Vulpiani, A , Phys. Rep., 156, 147 (1987). 

[5] R. Rammal, R., C. Tarmous, and A.-M.S.Tremblay, Phys. Rev. A, 31, 2662 (1985); R. Rammal, C. Tannous, P. Breton 
and A.-M.S.Tremblay, Phys. Rev. Lett. 54, 1718 (1985). 

[6] L. de Arcangelis, S. Redner and A. Coniglio, Phys. Rev. B, 31, 4725 (1985); Phys. Rev. B, 34, 4656 (1986). 

[7] A review for electrical noise can be found in: A.-M.S. Tremblay, B. Fourcade and P. Breton, Physica A, 157, 89 (1989); 
A.-M.S. Tremblay, and B. Fourcade, in Noise in Physical Systems, ed. CM. Van Vliet (World Scientific, Singapore, 1987). 

[8] A.W.W. Ludwig, and B. Duplantier, Phys. Rev. Lett. 66, 247 (1991). 



11 



[9] J. M. Deutsch, and R. A. Zacher, Phys.Rev. E 49, 8 (1994) and G.L. Eyink (to be published). 
[10] Y. Park, A.B. Harris, T.C. Lubensky, Phys. Rev. B 35, 5048 (1987). 

[11] B. Fourcade, PhD thesis, Universite de Sherbrooke, 1988. (unpublished). See also A. -M.S. Tremblay, and B. Fourcade, 
XVII th International Colloquium on Group Theoretical Methods in Physics, Eds. Y. Saint- Aubin and L. Vinet (Singapore, 
World Scientific, 1989). 

[12] The same structure also emerges in the context of dynamical systems , see : B. Fourcade et A.-M.S. Tremblay, J. Stat. 

Phys., 61, 607 (1990) and J. Stat. Phys., 61, 639 (1990). 
[13] A.-M.S. Tremblay, R.R. Tremblay, G. Albinet, and B. Fourcade, Physica A, 183, 398 (1992). 

[14] G. A. Garfunkel and M. B. Weissman, Phys. Rev. Lett. 55, 296 (1985). M. Octavio, J. Aponte and G. Gutierez, in Noise 

in Physical Systems, ed. CM. Van Vliet (World Scientific, Singapore, 1987). 
[15] B. Fourcade, P. Breton and A.-M.S. Tremblay, Physical Review B 36, 8925 (1987). 
[16] M.J. Stephen, Phys. Rev. B, 17, 4444 (1978). 

[17] R. Blumenfeld, Y. Meir, A. Aharony, and A.B. Harris, Phys. Rev. B 35, 3524 (1987). 
[18] A.B. Harris and T.C. Lubensky, Phys. Rev. B 35, 6964 (1987). 

[19] A. Erdelyi, W. Magnus, F. Oberhettinger, and F.G. Tricomi, Higher Transcendental Functions, Vol. II, p. 232-261, (McGraw 
Hill, New York, 1953). 

[20] G.N. Watson, A treatise on the theory of Bessel functions, p 306 (Cambridge, University Press, 1966). 
[21] G. Albinet, R.R. Tremblay and A.-M.S. Tremblay, J. Phys. I (France) 3, 323 (1993). 

[22] For a given order of the polynomial, say 2s, one can exploit the arbitrariness in the choice of homogeneous polynomials to 
select harmonic polynomials. These transform as irreducible representations of the rotation goup O(MN). See C. Itzykson 
and J.B. Zuber, Quantum Field Theory, p. 243 (Mc. Craw Hill, 1980). 

[23] B. Fourcade and A.-M.S. Tremblay Phys. Rev. B 39, 6819 (1989). 

[24] J.D. Brock, A. Aharony, R.J. Birgeneau, K.W. Evans-Lutterodt, J.D. Litster, P.M. Horn, G.B. Stephenson and A.R. 
Tajbakhsh, Phys. Rev. Lett. 57, 98 (1986). A. Aharony, R.J. Birgeneau, J.D. Brock and J.D. Lister, Phys. Rev Lett., 57, 
1012 (1986). 

[25] J.D. Brock, D.Y. Noh, B.R. McClain, J.D. Litster, R.J. Birgeneau, A. Aharony, P.M. Horn, and J.C. Liang, Z. Phys. B 
74, 197 (1989). 



12 



